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Abstract 

Using a 282 fb~^ data sample collected by the Belle experiment at the KEKB e^e^ collider, we 
study D^ decays to K~£~^i' and 7r~i~^i' final states. The D^ flavor and momentum are tagged 
through a full reconstruction of the recoiling charm meson and additional mesons from frag- 
mentation. The reconstruction method provides very good resolution in neutrino momentum 
and in q^ = (pe+Pu)'^- Normalizing to the total number of D^ tags, we measure the absolute 
branching fractions to be ^(L>° -^ Kiu) =(3.45 ± O.OTgtat ± 0.20syst)% and B{D^ -^ ttIu) = 
(0.255 lb 0.019g^at ^ 0-016gyg^)% and the semi-leptonic form factors (within the modified pole 
model) /f (0) = 0.695 ± 0.007gtat ± 0.022gygt and /^(O) = 0.624 ± 0.020gtat ± O.OSOgygt. 

PACS numbers: 13.20.Fc,14.40.Lb,13.66.Bc 



Exclusive semileptonic decays of B and D mesons are a favored means of determining 
the weak interaction couplings of quarks within the standard model because of their relative 
abundance and simplified theoretical treatment. The latter, given leptons are insensitive 
to the strong force, is due to the decoupling of the leptonic from the hadronic current. 
Limiting the precision on extractions of the couplings |Kib| and |V"cb| are our knowledge of 
the form factors parameterizing the hadronic current. Form factors from B and D meson 
semileptonic decay can and have been calculated usinglattice QCD techniques |lll2,|3| whilst 
heavy quark symmetry relates the two form factors \4^. Measurements of these decays are 
required to confront the theoretical predictions. In this Letter, we report measurements of 
the absolute rate and form factors of D° — *■ K~l~^i>i and D^ —>■ 7t~1~^i'i (l = e, fi), which have 
also been recently investigated by CLEO [s!, 'G^, BES l7| and FOCUS J8|. The measurement 
of D^ — >■ 7r~/i^z/^ is the first of its kind; furthermore, measurements of the form factor 
distributions f+{q^), where g^ is the invariant mass of the lepton pair, are substantially 
improved by using a novel reconstruction method with better g^ resolution than in previous 
experiments. 

Our analysis is based on data collected by the Belle detector |[9] at the asymmetric-energy 
KEKB storage rings jlOf with a center of mass (CM) energy of 10.58 GeV (T(4S')) and 60 
MeV below, corresponding to a total integrated luminosity of 282 fb^^. The Belle detector is 
a large-solid-angle magnetic spectrometer that consists of a silicon vertex detector, a 50-layer 
central drift chamber (CDC), an array of aerogel threshold Cherenkov counters (ACC), a 
barrel-like arrangement of time-of-fiight scintillation counters (TOF), and an electromagnetic 
calorimeter comprised of CsI(Tl) crystals (ECL) located inside a superconducting solenoid 
coil that provides a L5 T magnetic field. An iron flux return located outside of the coil is 
instrumented to detect K^ mesons and to identify muons. 

To achieve good resolution in the neutrino momentum and g^, we tag the D^ by fully 
reconstructing the remainder of the event. We seek events of the type e~^e~ — > DllgD*^X 
{D*~ -^ DgjgTT^}, where X may include additional n"^, n^, or K"^ mesons (inclusion of 
charge-conjugate states is implied throughout this report). Each candidate is assembled 
from a fully reconstructed "tag-side" charm meson (-Dtag) and additional particles {X), with 
the requirement that the combination be kinematically consistent with e^e^ -^ DllgD*^X. 
To the D^lgX is added a charged pion that is kinematically consistent with vr^ from D^^^ -^ 
DgjgVr". Candidate D[lgXiT~ combinations passing the analysis criteria thus provide a tag of 
Dgjg and its momentum without having detected any of its decay products. The decay D^^^ — *■ 
K^{7r^)i~V may thus be reconstructed with the neutrino momentum fully constrained. 

The Di*Jg is reconstructed in the modes D*+ -^ L)°7r+, D+7r° and D*° -^ D°7i°,D^j, 
with Z}+/° -^ i^^(n7r)++/+ {n = 1,2,3}. Each Z^tag and D^^^ candidate is subjected to a 
mass-constrained vertex fit to improve the momentum resolution. We require a successful 
fit of each Dtag candidate; furthermore, if this candidate is a daughter of a successfully fitted 
D*^g candidate, the event is treated as D^^^D*^X , otherwise it proceeds as Dtag-D*gX. The 
candidate X is formed from combinations of unassigned vr and K^K~ pairs, conserving 
total electric charge. The 4-momentum of D*" is found by energy-momentum conservation, 

assuming a Dt*gD*gX event. Its resolution is improved by subjecting it to a fit of the X 

tracks and the DIJ^ momentum, constrained to originate at the run-by-run average collision 
point, while the invariant mass is constrained to the nominal mass of a D*~ . The candidate 
is rejected if the confidence level of this fit is less than 0.1% (corresponding to =b3.3cr of 
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FIG. 1: Mass spectrum of D^- candidates. 



mass resolution). Candidates for 7r„ 



for each the candidate D^:„ 4-nionientuni is calculated from that of the Z)*„ and tt 



are selected from among the remaining tracks, and 

The 
corresponding invariant mass is shown in Fig. [1] The momentum is then adjusted by a 
kinematic fit constraining the candidate mass to that of the D°. For this fit, the decay vertex 
of the Dgj has been estimated by extrapolating from the collision point in the direction of the 
Dgj momentum assuming the average decay length; a comparison with the result without 
this vertex correction showed that the corresponding systematic error is negligible. Again, 
the candidate is rejected if the confidence level of this fit is less than 0.1%. 

Background lying under the D^jg mass peak (i.e. fake-Dgjg) is estimated using a wrong 
sign (WS) sample where the tag- and signal-side D candidates have the same flavor (-Dtag 
instead oi Dtpp-) ■ A MC study (including T(45') -^ BB and continuum (gg, where q = c, s, u, 
d) events [Ulll^) has found that this sample can properly model the shape of background 
except for a small contribution from real Dgj decays (~ 2%) from interchange between 
particles used for the tag due to particle misidentification. Background from fake D^^ 
is subtracted normalizing this shape in a sideband region 1.84 — 1.85 GeV/c^. We find 
56461 ±309stat±830syst signal D^jg tags. The systematic uncertainty derives from: statistics 
of the WS sample (0.5%); subtraction of real D^jg contamination in the WS sample (0.6%) 
and charge correlation in the background (2%). The latter was estimated with MC by 
comparing true right sign (RS) background with that of the WS. 

Within this sample of Dgjg tags, the semileptonic decay Dgjg -^ K^{Ti^)(i~V is recon- 
structed with K^{ti^) and l~ candidates from among the remaining tracks. An event is 
rejected if there are any remaining unassigned charged particles (2.2% [8.3%] of events in 
kaon [pion] mode) or if the remaining unassigned neutral energy exceeds 700 MeV (9.0% 
[12.0%] of events in kaon [pion] mode). This requirement has been optimized based on a 
comparison of the simulated and observed energy distributions; it removes a large fraction 
of D%„ hadronic decays with one of the final state charged tracks unreconstructed. The 
neutrino candidate 4-momentum is reconstructed by energy-momentum conservation, and 
its invariant mass squared, m^, is required to satisfy |m^| < 0.05 GeV^/c"^. 

Multiple candidates still remain in one third of Dg^g tags, and in about one quarter of 
the semileptonic sample. In these cases all candidates are saved and given weights adding 
up to 1. In the MC, the difference between the result of this method and the result when 
only the true signal is retained, is found to be negligible. 

The contribution from fake Dgj in the sample of semileptonic decay candidates is esti- 
mated using the Dgjg invariant mass WS shape of the D^jg tag sample, normalized in the 
previously defined sideband region. The effect of the additional selection criteria on the 



TABLE I: Yields in data, estimated backgrounds, extracted signal yields and branching fractions, 
where for the latter two, the first uncertainty is statistical and second is systematic; small differences 
in the numbers are due to rounding. 



channel 


full :D"ig 


K+e Ve 


K+pi v^ 


TT+e Ve 


TT+^^t'^ 


Yield 


95250 


1349 


1333 


152 


141 


fake D^ig 


38789 ± 830 


12.6 ±2.2 


12.2 ± 4.8 


12.3 ±2.2 


12.5 ±4.5 


semileptonic 


n/a 


6.7 ±2.6 


10.0 ± 2.5 


11.7 ±1.2 


12.6 ±1.9 


hadronic 


n/a 


11.9 ±5.6 


62.1 ±23.9 


1.8 ±0.7 


9.7 ±3.7 


signal 


56461 ± 309 ± 830 


1318 ±37 ±7 


1249 ±37 ±25 


126 ±12 ±3 


106 ±12 ±6 


Branching Fraction (10^'') 


345 ±10 ±19 


345 ±10 ±21 


27.9 ±2.7 ±1.6 


23.1 ±2.6 ±1.9 


(e and ^ channels, average) 


345 ± 7 ± 20 


25.5 ±1.9 ±1.6 



background shape has been conservatively estimated, by varying these criteria, to be 15% 
(35%) for e (/i) modes. 

There are also backgrounds from semileptonic decays with either an incorrectly identified 
meson or where additional mesons are lost in reconstruction. These backgrounds are highly 
suppressed by the good neutrino mass resolution. For D^- -^ 7r+£~z/ the most significant 
background is D^ -^ K^i^u amounting to 6% — 8% of the total yield. It was estimated 
using the reconstructed Z)° -^ K^i'u decays in data, reweighted with the (independently 
measured) probability of kaons to fake pions. Smaller backgrounds from D° -^ K*'^i~i' 
and D° — > p^i^u decays amounting to 0.8% — 0.9% were measured by normalizing MC 
to data in the upper sideband region ml > 0.3 GeV^/c^, which is dominated by these 
channels. For !D°g -^ K'^Crv, decays of U° -^ K*^Crv contribute at the level of 0.5%-0.8%, 
measured using a sideband evaluation as described above, while background from D^ — > 
TT "•"£"// and D^ -^ p^i^u was found to be negligible (< 0.07% of the total yield). Systematic 
uncertainties are assigned due to the following: MC statistics, which dominates overall (and 
according to channel ranges between 10% — 40% of the background); fake rate uncertainties 
(3% — 4%); and uncertainty on branching fractions of Dgjg -^ K^ /p'^i^u (^ 1%). 

The last source of background originates from D^- decays to hadrons, where a hadron 
is mis-identified as a lepton. It is measured with an opposite sign (OS) sample, where the 
lepton charge is opposite to that of the D*^ slow pion. Note that the signal is extracted 
from the same sign (SS) sample. In contrast to the SS sample, the OS sample has no 
signal or semileptonic backgrounds; fake D^jg are subtracted in the same manner described 
previously. Assigning well identified pion and kaon tracks a lepton mass, we construct pure 



background mf, distributions in both SS and OS, which we label /^ and /, 



fOS 



m 



K,7r. 



A fit of the weights ax and a^, of the components /^^ and f^^ in the ml distribution of 
the OS data sample is performed, and the hadronic background in the SS data sample is 
calculated as (ax/f-^ + cinfn^), utilizing the fact that the hadron misidentification rate does 
not depend on the charge correlation defining SS and OS. The method has been validated 
using MC samples. As the muon fake rate is about an order of magnitude larger than 
that for electrons, this background is much more significant for muon modes. Systematic 
uncertainties are assigned based on the bias of the method as studied in MC (11% — 35%), 
and parameter errors from the fit (16% — 35%). 



The signal yields and estimated backgrounds are summarized in the upper part of Table HI 
Efficiencies depend strongly on nx, defined as the number of 7r^*^°^ and K"^ mesons assigned 
to X (in e~^e~ — > dI^I,D*~X), and are determined with MC; typical ratios are ehiu/^no ~ 

& sig 

70%. As the observed nx distribution in the data is slightly different from that simulated, we 
reweighted the simulated efficiencies, amounting to a (+1.9 ± 3.9)% change in the efficiency 
correction; the corresponding uncertainty is accounted for in the systematic error. No other 
biases due to nx-dependent effects were found. 

Normalizing to the total number of D^- tags, the absolute branching fractions summa- 
rized in the lower part of Table H] are obtained; systematics are dominated by the absolute 
normalization. The results for electron decay modes are in good agreement with those 
from ja, lD], and the measured -D^jg -^ vr+e^z/ branching fraction confirms the pred iction 
of [l3|. The results for the muon modes are in agreement with the ratio given in 'l4l]. 

More information about the semileptonic decays is obtained by studying the differential 
decay width dV /dq^. The resolution in q^ of semileptonic decays is found to be aq2 = 
0.0145 ± 0.0007stat GeV^/c^ in MC signal events. This is much smaller than statistically 
reasonable bin widths, which have been chosen as 0.067 (0.3) GeV^/c^ for kaon (pion) 
modes, and hence no unfolding is necessary. Bias in the measurement of g^ that may arise 
due to events where the lepton and meson are interchanged, a double mis-assignment, was 
checked with candidate D^jg -^ K'^i~u events and found to be negligible. The differential 
decay width is bin-by-bin background subtracted and efficiency corrected, using the same 
methods described previously. 

In the theoretical description, the differential decay width is dominated by the form factor 
f+{<f) [l^- Up to order m| it is given by 

rirKM /^2 IT/ 12 



where Pk{tt) is the magnitude of the meson 3-momentum in the D^jg rest frame (lS|. These 
form factors have been calculated recently in unquenched lattice QCD [l|,|2]. In the modified 
pole model |15j], the form factor /_|_ is described as 

with the pole masses predicted as m{D*) = 2.11 GeV/c^ (for D^jg -^ K^i'u) and m{D*) = 
2.01 GeV/c^ (for D^jg — > Tc^i^u). Setting a^ = leads to the simple pole model jl6|. The 
ISGW2-model [H predicts 

where q^ax is the kinematical limit of q^ and aj is a parameter of the model, calculated for 
the i^-mode as ai{K) = 0.47 GeV-^c^. 

The measured q^ distribution is fitted with 2 free parameters to the predicted dif- 
ferential decay width dV/dq^ of the different models with /+(0) being one of the pa- 
rameters, and mpoie, Op or a/ respectively the other. Binning effects are accounted for 
by averaging the model functions within individual q^ bins. The fit to the simple pole 
model yields mpoie{K-i+v) = 1.82 ± 0.04stat ± 0.03syst GeV/c^ (xV^df = 34/28) and 



mpoie(vr-^+z/) = 1.97 ± O.OSstat ± 0.04syst GeV/c^ (xV^df = 6.2/10), in agreement with 
results from CLEO [5] and FOCUS J8|- While the pole mass for the 7r£z/ decay agrees within 
errors with the predicted value, m{D*), the more accurate fit oi mpoieiKiu) is several stan- 
dard deviations below m{D*). In the modified pole model, ap describes this deviation of 
the real poles from the m{Dt)) masses. Fixing these masses to their known experimental 
values, a fit of ap yields ap(D° -^ K-£+u) = 0.52 ± O.OSstat ± O.OGsyst (xV^df = 31/28) and 
ap(D° -^ 7T-i+u) = 0.10 ± 0.21stat ± O.lOsyst (xVndf = 6.4/10). Finally, a fit of the param- 
eter ai in the ISGW2 model yields a/(D° -^ i^-e+z/) = 0.51 ± O.OSstat ± O.OSsyst GeV'^c^ 
(xVndf = 33/28) and ai{D^ -^ -R-e+u) = 0.60 ± O.lOstat ± 0.09syst GeV-^c^ (xV^df = 
7.0/10) . Systematic uncertainties were studied using a toy MC where the exact simple 
pole model distributions for signal were randomly smeared according to the Gaussian errors 
found in the data. The fit reproduces the input pole masses without any significant bias; a 
shift of 1.2% (0.3(Tstat) observed in the pion mode was included in the systematic error. The 
subtracted background levels, which cause a correlation between g^ bins, were also varied in 
this toy MC. 

The fitted values for f\ '^{0) vary little for the different fits, for the modified pole model 
the results are f^{0) = 0.695 ±0.007stat ±0.022syst and /^(O) = 0.624 ± 0.020stat ±0.030syst; 
for the ratio (refitted without correlations due to normalization) we find 

f^fo)w\ ^ °'°^^ ^ 0.003stat ± 0.003syst (4) 

which is consistent within errors with the model-independent result using only the data in 
the first iriu q^ bin (g^ < 0.3 GeV^/c^). A recent theoretical prediction for the ratio fl] is 
0.040 ± 0.002 stat ± 0.005syst- Our resuh (JH) is in good agreement with those from CLEO ^ 
and FOCUS jlJ|, which measure slightly lower values. 

The measured form factors /+ ''^(q'^) are shown in Figure El with predictions of the simple 
pole model, unquenched j2| and quenched ^ LQCD. To obtain a continuous curve for /+ 
from the LQCD values reported at discrete q^ points, the values were fitted by a parabola, 
which is found to fit well within the stated theoretical errors and is not associated with 
any specific model. To quantify the degree of agreement, we calculate x^/^df between 
our measurement and the interpolated LQCD curve within the q^ range for which LQCD 
predictions are made. We find a x^/ndf of 28/18 (34/23) for the kaon modes and 9.8/5 
(3.4/5) for the pion modes; correlations induced by the fit of the calculated q^ points to a 
parabola have been considered. 

In conclusion, our measurement of the semileptonic 1^° -^ K{v:^lv decays yields absolute 
branching fractions in agreement with other new measurements, and the first measurement 
of the absolute branching fraction, B{D^ —>■ 7r~/i+z/). The good g^ resolution results in 
substantially improved measurements of D^ — * K^ {Ti~)e^ v /_|_(g^)-distributions. 
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FIG. 2: Form factors for (a) D° -^ R-i+u, in g^ bins of 0.067 GeVVc^ and (b) D^ -^ ir-i+u, in g^ 
bins of 0.3 GeV'^/c^. Overlaid are the predictions of the simple pole model using the physical pole 
mass (dashed), and a quenched (|3], light gray) and unquenched ([4], dark grey) LQCD calculation. 
Each LQCD curve is obtained by fitting a parabola to values calculated at specific q^ points. The 
shaded band reflects the theoretical uncertainty (without the BK-ansatz error for 2]) and is shown 
within the range of q'^ for which calculations are reported. 
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